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ABSTRACT 

Inherent binary or collective interactions in ensembles of quantum emitters induce a spread in the energy and lifetime of their 
eigenstates. While this typically causes fast decay and dephasing, in many cases certain special entangled collective states 
with minimal decay can be found, which possess ideal properties for spectroscopy, precision measurements or information 
storage. We show that for a specific choice of laser frequency, power and geometry or a suitable configuration of control 
fields one can efficiently prepare these states. We demonstrate this by studying preparation schemes for strongly subradiant 
entangled states of a chain of dipole-dipole coupled emitters. The prepared state fidelity and its entanglement depth is further 
improved via spatial excitation phase engineering or tailored magnetic fields. 


Introduction 

Ensembles of effective two-level quantum emitters consisting of single atoms, ions, or defects in solids are employed ubiq¬ 
uitously in quantum optics and quantum information. * They are the basis for precision spectroscopy or atomic clock setups, 
as well as for experiments testing fundamental concepts of quantum physics or implementations of the strong coupling cavity 
QED (quantum electrodynamics) regime.^’ ^ In the absence of direct particle-particle interactions, larger ensembles allow for 
faster, more precise measurements^ via a scaling of the effective single photon to matter coupling strength g by a factor y/N 
(with system size N) and a reduction of the quantum projection noise (by 1 / 

Eor any precise measurement one has to externally prepare, control and measure the particle dynamics. Hence, the emit¬ 
ters are almost unavoidably coupled to their environment. A suitable theoretical framework to model such experiments is 
open system dynamics with a coupling to a fluctuating thermal bath. At optical frequencies this can often be approximated 
by the zero effective temperature electromagnetic vacuum field. Still, extra perturbations by a thermal environment and 
background gas collisions cannot be avoided. 

In a laboratory experiment the particles need to be confined in a finite spatial volume that can be addressed by laser 
beams. Thus, increasing particle numbers will lead to higher densities, where direct particle-particle interactions as well as 
environmentally induced collective decoherence can no longer be neglected. Eor optical transition frequencies a critical density 
is conventionally assumed at the point where the average particle separation is of the order of an optical wavelength.® Above 
this limit vacuum fluctuations tend to become uncorrelated and decay becomes independent. However, recent calculations 
have shown that collective states can exhibit superradiance and subradiance even at much larger distances'** as long as the 
bandwidth of the emission is small enough. 

In many typical configurations and in optical lattices in particular, the particle-particle interaction is dominated by binary 
dipole-dipole couplings, with its real part inducing energy shifts and its imaginary part being responsible for collective de¬ 
cay."’*^ Generally, this interaction is associated with dephasing and decay. However, recently it has been found that under 
special conditions also the opposite can be the case and these interactions can lead to a synchronization or even a blockade 
of the decay. 

Oftentimes it is assumed that while such states exist, they cannot be prepared by lasers as they are strongly decoupled from the 
radiation fields. However, it was recently proposed that individual instead of overall addressing of the atoms can push the many 
particle system to evolve towards subspaces protected from decay or dephasing. When applied to Ramsey spectroscopy such 
states have been shown to exhibit frequency sensitivities superior even to those obtained from non-interacting ensembles.'® 
However, apart from special cases with an optimal lattice size and excitation angle, it is not so obvious how to implement such 
precise a control. 

In this work we highlight the surprising fact that interaction induced level shifts can be used to aid in preparing such states. 
In many cases the magnitude of the shifts a state experiences and its lifetime are tightly connected allowing one to identify 
and address interesting states via energy resolution. As a generic ensemble we particularize to a ID regular chain of quantum 
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Figure 1 . Selective state preparation procedure, a) A chain of N closely spaced quantum emitters (separation a with 
ka <^l, k being the laser wave number) are individually driven with a set of pumps {tj™}. b) The lasers are turned on for a 
time T, optimized such that an effective ;r-pulse into the desired subradiant target state is achieved, c) Level structure for the 
N systems where the -fold degeneracy of a given n-excitation manifold is lifted by the dipole-dipole interactions. The 
target states are then reached by energy resolution (adjusting the laser frequency) and symmetry (choosing the proper m). d) 
Scaling of the decay rates of energetically ordered collective states starting from the ground state (state index 1) up to the 
single- and double-excitation manifolds for 6 particles at a distance of a = 0.02 Aq. The arrows identify the decay rates for 
the lowest energy states in the single (A) and double (B) excitation manifolds, e) Numerical results of the time evolution of 
the target state population for A = 6 and a — 0.02 Aq during and after the excitation pulse. Near unity population is achieved 
for both example states A (where we used rj — 0.53 T) and B (for rj = 2.44r) followed by a subradiant evolution after the 
pulse time T shown in contrast to the independent decay with a rate T (dashed line). 


emitters coupled by dipole-dipole interactions with a tunable magnitude (by varying the interparticle separation). Collective 
coupling to the vacuum leads to the occurrence of subradiant as well as superradiant excitonic states.^** In particular, the 
subradiant states should prove extremely useful for quantum information as well as metrology applications as they exhibit 
robust, multipartite quantum correlations. As mentioned above, the atoms’ interactions provide a first handle for target state 
selection as they lead to energy resolved collective states. Furthermore, using a narrow bandwidth laser excitation matched to 
the target states both in energy and symmetry allows for a selective population transfer from the ground state via an effective 
Rabi ;r-pulse. 

In many cases, however, the required phase structure of the target state is not compatible with the excitation laser phase so 
that only a very weak coupling can be achieved. On the other hand, increasing the laser power reduces spectral selectivity 
by an unwanted addressing of off-resonant but strongly coupled states. Hence, to address a larger range of states of practical 
interest, we also propose and analytically study new methods of phase imprinting via a weak spatial magnetic field gradient. 
The small relative phase shifts increase the effective coupling to groups of emitters via a nonuniform phase distribution. With 
this method any state may acquire a finite laser coupling to the ground state via the magnetically induced level shifts resulting 
in an efficient population transfer with a minimal compromise on lifetime. 

The considered setup is a chain (see Fig. la) of N identical two-level systems (TLS) with levels |g) and \e) separated 
by a frequency of Oq (transition wavelength Aq) in a geometry defined by the position vectors {r,} for i = For 

each i, operations on the corresponding two-dimensional Hilbert space are written in terms of the Pauli matrices and 
raising/lowering operators < 7 ^ connected via af = < 7 A + 0 -r, = — i(( 7 A — 0 -r) and <7? = < 7 ^ 0 -r — < 7 r 0 -+. The complete 

Hamiltonian describing the coherent dynamics is 

H = Ho+Hdip = cooY,(^i^ + ^ 7 ’ ( 1 ) 

' A; 

where Hq is the free Hamiltonian and has degenerate energy levels (degeneracy = N\/{N — n)\n \ for level n) ranging from 
0 for the ground state to Aoo for the highest excited state. The second term Hdip describes interactions between pairs of TLS 
which can be induced either by an engineered bath (such as a common, fast evolving optical cavity field) or by the inherent 
electromagnetic vacuum. We denote the couplings between emitters i and j by and particularize to the case of a free-space 
one dimensional equidistant chain of TLS with small interparticle distances a such that a <C Ao (as depicted in Fig. la). 

For the sake of simplicity, we use dipole moments perpendicular to the chain for all numerical computations. To a good 
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approximation, in the limit of koa <C 1, the nearest-neighbor (NN) assumption can be used (such that Clij = Q.5ij±i) and exact 
solutions in the single-excitation manifold can be found. Within this subspace and approximation, the Hamiltonian assumes 
the form of a tridiagonal symmetric Toeplitz matrix with cOq on the diagonal and Q. above and below the diagonal. The solutions 
are readily available** with eigenvalues (Oq + e„, for an index m running from 1 to N, where e,„ = 20cos [nm/{N+ 1)] are the 
dipole-induced energy shifts. The corresponding eigenstates of the Hamiltonian are then 

I"*) = \G). with /- = y^l^sin , (2) 


where we used |G) = 

Spontaneous decay via a coupling to the free radiation modes in the evolution of the system can be included in a generalized 
Lindblad form,* 

-^[p] = ^ E Tij ). (3) 

^ ij 

where the Yij denote collective damping rates arising from the coupling to a common radiation field. These rates also strongly 
depend on the atomic distances a with two prominent limiting cases of Yij(a °o) = Fdij (independent emitters limit) and 
Yij{a —0) = r (the Dicke limit*^). In general, one can perform a transformation of the Liouvillian into a new basis by diag¬ 
onalizing the Yij matrix. This procedure leads to a decomposition into N independent decay channels with both superradiant 
(> T) and subradiant (robust) decay rates (< T).*® Note, however, that the states corresponding to these channels generally 
do not coincide with energy eigenstates of the Hamiltonian, so that we cannot reduce the system dynamics to simple rate 
equations. 

Results 

Selective state preparation 

Tailored coherent excitation. As mentioned above, our dipole coupled systems possess states with a large range of radiative 
lifetimes and energy shifts. Depending on the desired application particular states can be highly preferable over others. In a 
first straightforward approach we now illustrate that in principle it is possible to access a desired collective state simply by a 
selective coherent driving with a properly chosen amplitude and phase for each TLS. This is described by the Hamiltonian 

Hm = '£ ( 4 ) 

j 

with a suitably chosen set of VjJ. For a targeted eigenstate in the single-excitation manifold, some analytical insight on how 
to choose these amplitudes can be gathered from the state’s symmetry. For energy eigenstates this can be found quite reliably 
within the NN approximation.^** In an equidistant finite chain our calculation suggests the following choice of driving fields 
at laser frequency O)/, 


77 ™ = 77 sin 


/ nmj \ 

U+iy’ 


(5) 


chosen to fit the symmetry of a target state \m). 

The selectivity of the excitation process can be further improved by an energetically resolved excitation of a given state 
\m) by a proper choice of the laser frequency (Oi = (O 0 + £m its bandwidth. This is possible due to the interaction induced 
level splitting from H^np (as depicted in Fig. Ic). Indeed, in perturbation theory and in a frame rotating at (Oi the evolution of 
the system starting from the ground state up to a normalization factor leads to 

~ |G)-7T7f|m). (6) 

The success of the corresponding process is illustrated in the sequence of plots in Fig. 1, where the \m = N) state with n—1 
is considered (target state A) and accessed via the combination 77 ^ of pumps lasting for a duration T. 

Numerical simulations were performed on a six-atom chain with driving strength 77 = 0.53 F at an interatomic separation 
of a = 0.02 Aq. The time for which the pumps are switched on is T = l.SSF^* which is considerably shorter than the time 
scale governed by the decay rate of 0.0009 F of the target state. The resulting dynamic is an effective ;r-pulse (efficiency of 
99.94%) flipping the population into the state \m = N) followed by an extremely slow decay, indicating the robustness of the 
target state (as seen in curve A of Fig. le). 
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It is, of course, desirable to target higher excitation manifolds as well. In the absence of analytical expressions or good 
approximations for the target states, we employ phases that yield maximal asymmetry, i.e. fjj = T 7 (—1)-' for any j = 1,...,A^. 
Such a driving can be expected to address collective states, where the helds emitted by any two neighboring particles interfere 
destructively^'* (similar to a previously investigated mechanism*^). Numerical simulations show that the resulting collective 
states indeed exhibit the lowest energy shifts of the targeted manifold and can be expected to be long lived. The resonance 
condition for a specihc state jy/) within the manifold n is ncoi = ncoo + dco^,, where 5(0^ = {\j/\Hdip |y/). As an illustration, 
the curve B in Fig. le shows an almost perfect efficiency (98.36%) two-photon ;r-pulse allowing for a population transfer to 
the longest-lived collective state in the second excitation manifold of A = 6 emitters separated by a = 0.02 Aq. The chain was 
driven with a strength of t] = 2.44T for a time T = 3.44r^*, which again is significantly shorter than the natural time scale 
given by the target state decay rate of 0.0402 T. 

Let us add a comment on the practical implementation of such an addressing. In typical current experimental configura¬ 
tions for clocks based on ID magic wavelength lattices^*’^^ the atoms are very close and hardly allow for an individual direct 
particle addressing. One is largely limited by a quasi plane wave driving, which typically addresses all particles with equal 
intensity. If the pump light is applied perpendicularly to the trap, the evolution is governed by a symmetric Hamiltonian Hsym, 
obtained from equation (4) with an equal pump amplitude T]J = r] for any m and j. A laser excitation from the ground state 
into the state \m) is connected to the coupling amplitude Xm = {fn\Hsym \ G) = V which yields 

{ 0 if m is even, 

ifis odd. 

We will refer to states with even m as dark states as they cannot be accessed by the laser excitation and call the remaining 
ones bright states. In the limit of large atom numbers A ^ 1, it is of interest to investigate the two cases, where m <^N 
and m ^ A, for states at the top/bottom of the manifold. In the hrst case, the function for the driving yields Xm ~ t? v/8A /mn, 
whereas in the other case we have Xm ~ 0. 

Note, that sometimes geometry can change this behavior. For a ID string of equidistant emitters illumination at a chosen 
angle of incidence and polarization leads to a designable phase gradient of the excitation amplitudes. The situation becomes 
even more complex for a 3D cubic lattice, where the phases also differ in the different lattice planes. As a lucky coincidence, a 
perpendicular plane illumination at the clock frequency in a magic lattice for Strontium (Sr) targets an almost dark state. This 
leads to subradiance and in principle allows for a spectral resolution better than the natural linewidth.^^ In not so favorable 
cases one could also think of a specific lattice design to facilitate a tailored dark state excitation. 

Radiative properties. In order to be useful resources for quantum information applications, target states should exhibit 
robustness with respect to the environmental decoherence. To identify states of minimum decay rate, we scan through the 
eigenstates | V4) of the Hamiltonian H ~ Ho+ Hdip (for k= \ and compute their decay rates F (see section Methods 

below). We find that generally, for a given manifold, the energetic ranking of the states closely indicates their robustness to 
decay (as illustrated by the color-coding in Fig. Ic) ranging from blue for subradiant states to red for superradiant states. This 
is due to the fact that both radiation and energetic shifts are strongly dependent on the symmetry of the states. In Fig. Id, for 
A = 6, we plot the decay rates of the collective states in the hrst (n = 1) and second (n = 2) excitation manifold arranged as a 
function of their increasing energy corresponding to the level structure of Fig. Ic. Superradiant states are found at the upper 
sides of the manifolds while the ideal robust states lie at the bottom. In Fig. Id, the arrows indicate the optimal decay rates 
in the single- (0.0009F) and double-excitation manifolds (0.0402F) corresponding to target states A and B whose population 
evolution is depicted in Fig. le. 

Within the single-excitation manifold, an analytical expression for the decay rate of a state \m) can be found as Ym = 
T.i.jYijfrfj'- For small distances the state m = 1 (upper state) is superradiant, whereas states at the bottom of the mani¬ 
fold N exhibit subradiant properties. In the Dicke limit where a = 0 we have jij = F for any i and j, and we can compute 
r„, = 2r cot^ [mn/ (2A-|-2)]/(A-f 1) for m odd and r„, = 0 for m even. Note, that in this particular limit, these are the same 
conditions as for the darkness and brightness of a state. For large numbers of emitters, we recover the expected superradiant 
scaling with A for the state with m = 1, i.e. Fi Ri SFA/ n^. On the other hand, large m yield a decay rate of Y,„ k, 0 (perfect 
subradiance) in the same limit. 

There are two important conclusions from these results: i) since in the considered limit the decay rate of the superradiant 
state \m = 1) scales with Fi A, whereas its driving is X\ Vn, driving this state becomes more difficult with increasing 
atom number due to the reduced time-scale and ii) if the number of atoms is not too large, Xm will remain finite, while r„, 
already indicates vast subradiance due to its scaling-down with A. Hence, there are robust states that remain bright, i.e. they 
can be driven directly even though the driving is not matched to their symmetry. 
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Accessing dark states via magnetic fieid gradients 

The direct symmetric driving with Hsym allows access to bright states only. Given that nearby dark states can conceivably be 
more robust, we now employ a progressive level shifting mechanism that allows for a coupling between bright and dark states. 
This is achieved by subjecting the ensemble to a magnetic field with a positive spatial gradient along the chain’s direction. The 
increasing energy shift of the upper atomic levels (as depicted in Fig. 2a) plays a role similar to the individual phase imprinting 
mechanism described previously. For each particle the shift of the excited level induces a time-dependent phase proportional 
to the value of the magnetic field at its position. We demonstrate the mechanism for a particular two-atom example, where 
indirect near unity access to the dark subradiant asymmetric collective state is proven and extend it to the single-excitation 
manifold of N atoms. 

Two-atom case. The eigenstates of the Hamiltonian 
Ho+Hdip are jF) = \ee), |G) = |gg) and in the single¬ 
excitation subspace |5) = {\eg) + \ge))/^/2 and |A) = 

(k^) ~ k^))/'\/2- The symmetric state |5) is superradi¬ 
ant (F^ = Fi = F -f 712 ) and bright, directly accessible via 
symmetric driving with strength X\ = v/2t?- The asym¬ 
metric state |A), on the other hand, is subradiant (F^i = 
r 2 = F — 712 ) and dark. Indirect access can be achieved 
by shifting the second atom’s excited state by 2 Ab (see 
schematics in Fig. 2b), where Ag is tunable and quantifies 
the per-emitter shift for a given magnetic field amplitude. 

We first analyze the dynamics in the absence of decay by 
solving the time-dependent Schrodinger equation governed 
by the Hamiltonian H = Hq + H^ip -f Hsym + Hb, where 
Hb = IAbO^ O 2 ■ We reduce the dynamics to three states, 
and assume a quasi-resonant Raman-like scheme where the 
population of jF) is at all times negligible. An effective 
two-level system arises (between the ground state and the 
asymmetric state; see section Methods below) and the res¬ 
onance condition can be identified as 




2n 


■|A) 


.|G) 



tP-') = -Ab + \/aI -7 - 2tj2, 


( 8 ) 




rt 


with an effective Rabi frequency of 

V® =-(9) 

n-t-Y^Al -7 02-2772 

To fulfill jc^k <c 1, we need to restrict the driving to a pa¬ 
rameter regime where T 7 ,Ab ^ O. A scan over the mag¬ 
netic field is performed and the exact numerical results 
for the asymmetric state population are plotted in Fig. 2d 
against the adiabatic solution showing near unity popula¬ 
tion transfer for an optimized A^. Further restrictions are 
imposed when decay is considered. These stem from the 
fact that the coherent process described by Vr should be 
faster than the incoherent one characterized by F^. For 
close particles, the ability to tune the distance ensures 
that the scaling down of r ,4 is very fast and the above 
conditions are readily fulfilled. For the particular exam¬ 
ple illustrated in Fig. 2d we chose a — 0.05 Aq, resulting 
in fi = 23.08r, r ,4 = O.OIOF. The 0.994 population is 
reached atT = 16.19r^*, which is very close to the theoret- 
ical estimate of T = n/2Vp'^ = 16.179r^^ obtained from 
the adiabatic solution under the assumption of a ;r-pulse 
transferring the population to the target state. 


Figur© 2. Coupling to dark states via a magnetic field 
gradient, a) Linearly increasing level shifts along the chain 
occuring in the presence of the magnetic field gradient, b) 
Illustration of the level structure and indirect dark state access 
for two coupled emitters. While symmetry selects the state \S), 
off-resonant addressing combined with bright-dark state 
coupling of strength Ar allows for a near-unity population 
transfer into the state |A). c) Dynamics in the single-excitation 
manifold of N coupled emitters where symmetric driving 
reaches the bright states with amplitudes Xm while the 
magnetic field couples neighboring dark and bright states, d) 
Plot of the asymmetric state population for the two-atom case 
as a function of the increasing magnetic field (solid line) 
compared to the steady-state approximation (dashed line) at 
numerically optimized time T = 16.19r^*, with parameters 
7] = F and a — 0.05 Aq- e) For a chain of = 4 emitters, a 
91 %-efficient tt- pulse to the most robust state can be achieved 
as demonstrated in the population evolution plot. The 
separation is chosen to be a = 0.025 Aq, while 77 = 40r and 
numerical optimization is employed to find Ab = 0.98 F. 
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Many-atom case. For a chain of N atoms, we consider the 

progressive shifting of excited levels along the chain depicted in Fig. 2a. This is realized by the application of a magnetic field 
with a constant gradient and is described by the Hamiltonian Hb — 1)(7^(7;^. Let us consider a dark state \d) {d 

even) and the bright state \b = d—1) immediately above. Their coupling via Hb is quantified by A^b = 2AB^,(i— l)/J^/f, as 
shown in Fig. 2c. 

We develop a protocol where direct off-resonant driving into the bright state (amplitude Xb) combined with a coupling between 
the bright and dark states via the magnetic field leads to an almost unity population transfer into the dark state. Given a 
sufficient energy separation, the problem can be reduced to solving the time-dependent Schrodinger equation for the three 
coupled state amplitudes Cb,Cd and eg. Following the same adiabatic approximation as in the two-atom case we reduce the 
general dynamics to an effective two-level system between the states meant to be connected by an effective TT-pulse, i.e. \d) 
and |G). The generalized resonance condition (with e^b = £d — £h) reads 


At^) 


-AB(iV-l) 


£d + £b 



/y 2 


( 10 ) 


and was obtained in the limit where the coupling of the dark state to the other adjacent bright state \d+ 1) was neglected 
owing to the relation Xd-i ^ Xd+\- The effective transition rate between the ground state and the state \d) is 


(iV) ^ Xb\^db\ 

^ A -f £/, -f Ab {N — 1) 


( 11 ) 


The addition of decay imposes a new constraint on the timescale of the process, i.e. ^ Fj, required to ensure near unity 
population in the dark state. The fulfillment of this condition depends on the individual system under consideration. As an 
illustration of the procedure. Fig. 2e presents the targeting of a robust dark state in the single excitation manifold of four 
particles. Note, that the numerical results are performed in an exact regime beyond the NN approximation and are in excellent 
agreement with our conclusions obtained from the NN treatment. 


Discussions 
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Entanglement properties 

To justify the usefulness of collective states for quantum information purposes, we employ the von Neumann entropy to 
analyze their entanglement properties. More specifically, we compute the von Neumann entropy of the reduced density 
matrix of a single two-level emitter (showing the degree of its bipartite entanglement with the rest of the system) de¬ 
fined by S{ps) = — where A, is the i-th eigenvalue of p^ and 01og20 = 0. We furthermore minimize the set of 

values for all atoms to obtain a lower 
bound on the entanglement contained 
in the system. We compare the nu¬ 
merical results to the single-atom en¬ 
tropy of the symmetric Dicke state 
\—N/2,—N/2 + n)}^ For these partic¬ 
ular states the entropy is maximized if 
the number of excitations in the state is 
n = N/2. It follows that it is highly de¬ 
sirable to drive the system into robust 
states as close as possible to n = \N/2\ 
excitations (where \N/2\ is the largest 
integer smaller or equal to N/2), since 
this manifold contains the most entan¬ 
gled state. A comparison of the exact 
numerical data and the analytical expres¬ 
sion for the entropy is shown in Fig. 3a. 

Another way to characterize the 
entanglement of the prepared state is 
to investigate their depth of entangle¬ 
ment ,which does not quantify the 
entanglement itself but rather shows 
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Figure 3. Entanglement properties, a) Comparison of the numerically computed 
von Neumann entropy (empty circles) of the reduced density matrix of the chain 
minimized over the atom index and the analytical expression for the entropy of 
the Dicke state (green circles), both for excitations n— I and n = \N/2\ as a 
function of the atom number N at distance a — 0.1 Aq. b) Depth of entanglement 
of the subradiant four-atom state (blue dot) prepared by the magnetic field 
gradient scheme (see Fig. 2e). It clearly lies above the k = 2> boundary indicating 
four-atom entanglement. The k-atom entanglement boundaries of the target state 
population Pt as a function of the ground state population Pq have been computed 
for the corresponding target state of a four-atom chain at distance a = 0.025 ^ 













how many atoms of an ensemble are 
involved in the present entanglement. 

This measure has been used in recent ex- 

periments^^’^® since it is a readily measurable quantity. The depth of entanglement is computed as follows: given an A^-atom 
target state in which an arbitrary number of said N atoms is entangled, we compute the limit of how much population one can 
drive into this state such that the resulting density matrix p remains separable into a subset of density matrices that exhibit 
no more than A:-atom entanglement (1 < A: < N). This may be done by numerically maximizing the target state population P, 
as a function of the ground state population Fq for different k. The boundaries themselves indicate how many atoms need to 
be entangled in order to prepare the pure target state, i.e. the boundary where the target state population is maximized to 1 
corresponds to the number of atoms entangled in the (pure) target state. If a general prepared state has a target and ground 
state population such that the corresponding data point lies on or above the A:-atom boundary, more than k atoms are entangled. 

Obviously, for the pure target states considered in the above computation all atoms contribute to the entanglement, since 
otherwise the minimal von Neumann entropy as shown in Fig. 3a would be zero. For a more interesting result, we can compute 
the depth of entanglement in order to demonstrate the efficiency of the driving procedure using a magnetic field gradient as 
in Fig. 2e. From Fig. 3b, where all boundaries have been plotted for the considered subradiant four-atom state, it is clear that 
the prepared state shows all-atom entanglement as the corresponding data point lies far above the boundary for three-atom 
entanglement. 

Implementation considerations 

The proof-of-principle technique presented above has been particularized on a specific generic system of emitters in an equidis¬ 
tant chain. The choice is natural since the electromagnetic vacuum provides a simple example for both collective dispersive 
and dissipative dynamics. To exemplify a possible realization we consider a particular system^^ where bosonic Sr atoms are 
trapped in a magic wavelength optical lattice at separations of a = 206.4nm. The working transition is at Aq = 2.6/rm, between 
the ^Pq and ^Di electronic states. This amounts to a ratio offl/Ao~ 1/13 which allows for an operation in the regime targeted 
by our scheme. The corresponding single atom decay rate is at the order of F = 0.3 MFIz and circularly polarized light can 
allow for transitions between states with a difference of 1 in magnetic quantum number. We have numerically investigated a 
system of 4 atoms in such a configuration and found a sizeable 73% target state population for rj = IF and Ab — 0.5 F, under 
the conditions of a relatively small level shift between the dark and bright state around 6 F which does not allow for large 
driving powers. For further optimization of the efficiency of the target state preparation one could envision a modified setup 
where a trapping transition of smaller wavelength can be chosen that would most importantly allow for better state separation 
(owing to larger dipole shifts). The corresponding magnetic field gradient required to produce the considerable = 0.5 F 
shift on a distance of a = 206 nm is around 5.2 • 10^ G/m, not far from state-of-the-art values achievable in high magnetic 
field gradient magneto-optical trap experiments.^^’^® Of course, there are many detrimental practical effects that can seriously 
limit the above technique such as light-assisted collision loss. We envision the extension of the described technique to systems 
where both the coherent and dissipative particle-particle interactions can be suitably tailored. For example, the same kind 
of dipole-dipole Hamiltonians can occur in 3D lattices of polar molecules^*' or between two different color NV centers in 
diamonds.^* 

Conclusions 

Direct particle interactions are typically detrimental and limiting in precision measurement applications. Here, we have 
presented some specific opposite examples, where the collective nature of the decoherence combined with the coherent binary 
dipole-dipole interactions is used as a new resource for the controlled and efficient preparation of specially selected states. The 
excitation scheme can be tailored to address target states exhibiting both entanglement as well as robustness against decay. As 
a generic example we studied the case of a one-dimensional system of tightly spaced equidistant quantum emitters. Already 
the inherent dipole-dipole coupling allows for a targeted state preparation technique via energy selection. The performance 
of the excitation can be enhanced additionally via the continuous application of a spatially increasing magnetic field. The 
general principle of such a phase imprinting technique is potentially applicable in many specific environments such as optical 
lattices or atoms and ions localized within one or more common optical cavity modes,^^’^^ NV-centers or superconducting 
qubits coupled to CPW transmission lines or resonators. 

Methods 

Decay rate of the states 

In order to arrive at an analytical expression for the decay rate of an eigenstate |v 4 ) of the Hamiltonian in equation (1), we 
consider the homogeneous part of the differential equation of the corresponding density matrix element that arises from the 
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master equation. The solution of this differential equation yields an exponential decay. The rate at which the state population 
decays may be written as 

^¥k = -{y'k\-^[Wk){¥k\]\¥k)=Y.yu{y'k\(yi^(yJ\yfk)- ( 12 ) 

Note, that this is true only for states that contain one specific number of excitations, i.e. they are eigenstates of the operator 
£,■<7?. Obviously, this is fulfilled for eigenstates of the considered Hamiltonian. Equation (12) was used in order to compute 
the rates depicted in Fig. Id and throughout the manuscript. For example, we used it in order to compute the decay rate of the 
eigenstates in the NN approximation Fm- 

Subradiance and disorder 

Fet us consider the influence of positioning disorder on subradiant properties of the target states. To mimic disorder we perturb 
an equidistant chain of N emitters (average separation a) by introducing an uncertainty in each emitter position quantified by a 
defect parameter s (normal distribution of variance sa). We then write the randomized matrix of decay rates and find the mini¬ 
mum decay channel without as well as in the presence of disorder of i = 20% and s — 40%. For the s = 0% case, it has been 
shown^® that the minimum decay rate scales exponentially with N even for distances up to 0.4Ao, while the linear scaling with 
N typical for superradiance is reached for a <C Ao only. After averaging over 100 random configurations, we plot the logarithm 
of the minimal rates 
As a somewhat surprising result, subra¬ 
diance scales even better with N as the 
disorder increases. This might be under¬ 
stood as a destructive interference effect 
brought on by the cancelation of emit¬ 
ted photons stemming from the random 
positioning. As pointed out in previous 
investigations,^® the states of low sym¬ 
metry (as, for example, the m=N state) 
possess decay rates closest to the analyt¬ 
ically derived minimal rate. We analyze 
the respective sensitivity of the state sub¬ 
radiance to disorder by initializing the 
system of N emitters in the m—N state 
and allow it to decay. The outcome is 
plotted in Fig. 4b and shows remarkable 
robustness of the disordered systems on 
a long time-scale. While on a short 
time-scale disorder pushes the consid¬ 
ered state into faster decaying channels, 
the long time limit shows that the re¬ 
maining population accumulates in the 
disorder-enhanced robust states. 

For short time-scales, the state still decays slowly (subradiantly), however, the decay rate increases with growing disorder 
(s = 40%). More remarkable, though, is the behavior the decaying states show for long time-scales, as the states subject to 
larger disorder become more robust than the unperturbed system. This is due to the fact that all population in the m = N 
state that decays through more radiative channels have decayed at that point and only the most subradiant channel (minimal 
eigenvalue of the decay rate matrix) remains. As seen in Fig. 4a, this eigenvalue is even further reduced by disorder which 
explains the long time-scale behavior in Fig. 4b. 

Coherent dynamics with a magnetic fieid gradient 

Two-atom case. To find the expressions in equation (8) and equation (9) we solve three coupled differential equations neglect¬ 
ing the population of the fully inverted state \E) as far off-resonant for all times. In the collective basis, where any state may 
then be written as 11^) = |5) -f ca |A) +cg\G), the equations are 


ids = {A + Ab + Fl)cs - Abca + V2r]CG, 

(13) 

icA = {A + Ab- Q)ca - Abcs, 

(14) 

icG ^rjcs, 

(15) 


as a 

function 

of 

increasing 

N 

in 

Fig. 

4a. 

a 



b 








Figure 4. Subradiance and disorder, a) Plot of the logarithm of the minimal 
eigenvalue of the decay rate matrix (matrix with entries jij) as a function of N at 
a distance of a = 0.4Ao for increasing levels of disorder (s = 0,0.2,0.4). b) 
Decay of the \m = N) state as a function of time. In the presence of disorder 
(s = 0.2,0.4) the short time and long time behaviors are fundamentally different. 
At short times, disorder can push the state towards faster decaying channels 
while decay inhibition due to disorder occurs at larger times. 
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where Q. = i 2 i 2 is the coherent interaction between the atoms and A is the detuning between the atomic resonance frequency 
and the driving laser. For an efficient driving of |A) the population of the state \S) needs to be negligible which allows us to 
set a steady-state condition, namely C 5 = 0 yielding the desired effective two-level system between |G) and |A). 

Many-atom case. The same approach as in the two-atom case may be used to describe the dynamics in the single-excitation 
manifold for an arbitrary number of atoms in a chain. Given sufficient energy separation we may neglect all states but the 
ones we aim to address. We can indirectly address a dark state \d) by driving the bright state \b) immediately above, which 
is coupled to the dark state by a magnetic field gradient. Neglecting all populations but c&, Cd, and cg and their respective 
couplings via the magnetic field gradient, the investigation reduces to the equations 


ich — [A-|-e;,-|-AB(N— i)]ch + AdbCd + XbCg, 

(16) 

icd = [A -f -f Ab(N— 1)]c^ + /^dbCb^ 

(17) 

icG = XbCb- 

(18) 


For an efficient driving of the dark state we may again invoke a steady-state condition on the bright state population Cj, = 0. 
This, again, yields an effective two-level system between the ground and the dark state with resonance condition and Rabi 
frequency as displayed in equation ( 10 ) and equation ( 11 ), respectively. 

Von Neumann entropy 

For a Dicke state an analytical expression for the von Neumann entropy of the reduced density matrix can be obtained. First, 
note that, since Dicke states are invariant under a permutation of the atoms, all reduced density matrices are identical. Hence, 
they all share the same von Neumann entropy for a given number of excitations n. We may choose to reduce the full density 
operator p to the density matrix of the first atom in the ensemble, i.e. p} = Ps = ..,n{p) which yields a von Neumann 

entropy of 

S(ft) = ^logj (^) - (1 -^) log2 (1 - (19) 

For the actual eigenstates of the Hamiltonian in equation (1) this computation needs to be done numerically. Furthermore, 
these states are not invariant under permutation of atoms and hence it is required to minimize the entropy with respect to the 
atomic chain index in order to find the lower bound. 


Depth of entanglement 

The boundaries depicted in Fig. 3b were found by maximizing the target state population with the condition on the density 
matrix of the prepared state to contain no more than k-atom entanglement, i.e. p — p,' with k, < k and at least one k, = k. 
To compute the boundaries we generalized the algorithm that was previously used solely for the W-state^^ to arbitrary states 
in the single-excitation manifold. For the computation of all boundaries we need to distinguish the two cases where Pg = 0 
and Pg > 0. Considering a separable state (k = 1), the boundary for Pg>Q is found to be 


max{Pt) = Pg max 
n, ai=Vfb 


Ll‘^'1 

i 


Oi 


2 


( 20 ) 


where a,- S [0,1] and c, are the coefficients of the target state. For Pg = 0 the maximization is much simpler, i.e. max(P,) = 
max I c, I , which is found by setting one a, = 0 and the remaining coefficients = 1. Note, that for both these and all 
following computations we neglect the symmetry of the state, i.e. the phases of the coefficients a by using |c,|. This is valid 
due to the invariance of entanglement under local unitary operations and necessary if we restrict the coefficients a, in the way 
we did. 

For multiple-atom entanglement (k > 1) the matter of finding the corresponding boundary is no longer so simple. In order 
to find the maximum population, we assume maximally allowed entanglement in the prepared state. We split the prepared 
state into M = \N/k\ sets, where M—\ sets are k-atom entangled and the remaining one isk! = N — k{M — 1 )-atom entangled. 
To find the maximum, one has to consider all possible positions of the k'-entangled state. If, for example, the k'-entangled 
state is at the last position, the population of the target state |f) in the prepared state reads 


P,= 


{t\ 


/M-l 


l<pf) ® I<Pm) 


i=\ 


( 21 ) 
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where 


|(pf) = Ui \Gk) + 


1 - L \^k) 


( 22 ) 


is a general non-separable state of k atoms in the single-excitation manifold. The state |G^) is the A:-atom ground state and the 
coefficients G [0,1] have to be normalized, i.e. = 1 V i. One then has to maximize the target state population with 

respect to the coefficients a, S [0,1] and Xj. with the condition H; OCj = The number of these coefficients, however, 

grows vastly with the number of atoms, hence numerical computations are limited. For Pq — 0 one can again choose one 
a,- = 0 and all , = 1. 

Note, that all boundaries computed via this maximization only hold for pure states. In order to find the boundaries for 
mixed states we need to compute the convex hulls of the respective boundaries.The k = N boundary is found when a perfect 
superposition between the ground and target state is reached. 

In this work we considered the specific case of an exciton state of a four-atom chain. In that case, when investigating 
two-atom entanglement the permutation of the k'-entangled state is rendered unnecessary since k' — k = 2. Unfortunately, this 
is no longer true for k —3, where we did have to account for all permutations. 
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